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Abstract 

Background: Nucleosomes are present throughout the genome and must be dynamically regulated to 
accommodate binding of transcription factors and RNA polymerase machineries by various mechanisms. Despite 
the development of protocols and techniques that have enabled us to map nucleosome occupancy genome-wide, 
the dynamic properties of nucleosomes remain poorly understood, particularly in mammalian cells. The histone 
variant H3.3 is incorporated into chromatin independently of DNA replication and requires displacement of existing 
nucleosomes for its deposition. Here, we measure H3.3 turnover at high resolution in the mammalian genome in 
order to present a genome-wide characterization of replication-independent H3.3-nucleosome dynamics. 

Results: We developed a system to study the DNA replication-independent turnover of nucleosomes containing the 
histone variant H3.3 in mammalian cells. By measuring the genome-wide incorporation of H3.3 at different time points 
following epitope-tagged H3.3 expression, we find three categories of H3.3-nucleosome turnover in vivo: rapid turnover, 
intermediate turnover and, specifically at telomeres, slow turnover. Our data indicate that H3.3-containing nucleosomes 
at enhancers and promoters undergo rapid turnover that is associated with active histone modification marks including 
H3K4me1, H3K4me3, H3K9ac, H3K27ac and the histone variant H2A.Z. The rate of turnover is negatively correlated with 
H3K27me3 at regulatory regions and with H3K36me3 at gene bodies. 

Conclusions: We have established a reliable approach to measure turnover rates of H3.3-containing nucleosomes on a 
genome-wide level in mammalian cells. Our results suggest that distinct mechanisms control the dynamics of H3.3 
incorporation at functionally different genomic regions. 



Background 

Incorporation of histone variants into chromatin critic- 
ally influences the properties of nucleosomes that play 
important roles in regulating transcription and epigenetic 
memories [1-6]. The histone variant H3.3 differs from 
canonical H3 by a few amino acids and is ubiquitously 
expressed in eukaryotes. Different from canonical histone 
H3, which is expressed in S phase and is incorporated into 
chromatin during DNA replication, H3.3 can be incorpo- 
rated into chromatin independent of DNA replication 
[7,8]. The incorporation of histone variants is tightly regu- 
lated by histone chaperones. The H3.3-specific chaperones 
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Atrx-Daxx and HIRA deposit H3.3 primarily at telomeres 
and non-heterochromatic regions, respectively [9-12]. Des- 
pite small differences in amino acid sequence between 
H3.1, H3.2 and H3.3, each is distributed differentially across 
the genome and carries its own characteristic histone modi- 
fication signature, strongly suggesting distinct functional 
roles for each H3 variant [9]. Histone marks that are 
associated with gene activation, such as acetylation marks 
and H3K4me3, are typically found on H3.3 whereas 
H3K27me2 and H3K9me3 are found on H3.2 [13]. Marks 
associated with gene silencing are predominantly found 
on H3.1. However, the precise relationship between H3.3 
deposition and transcription are not well understood. In 
Drosophila, H3.3 replaces H3 during gene activation and 
becomes deposited in active chromatin and especially very 
highly expressed ribosomal DNA [14,15]. In mammalian 
tissues the pattern of H3.3 enrichment is also associated 
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with gene activity and H3.3 is generally associated with 
the transcription start site (TSS), transcription end site 
(TES) and gene bodies of active genes [16], although the 
unique chromatin of embryonic stem cells (ESCs) also 
carries H3.3 at promoters of certain inactive genes [1,9]. 
In spite of their independent evolution, plant H3.3 and 
H3.1 display a broadly similar distribution to animal H3 
variants, which indicates a conserved function for H3 
variants [17]. Nucleosome occupancy and positioning are 
critical to regulating gene transcription and epigenome 
maintenance [18,19]. In mammalian cells, although nu- 
cleosomes are present at promoters and enhancers, they 
must be dynamically regulated to accommodate binding 
of transcription factors and RNA polymerase machineries 
by various mechanisms [20-23]. 

Intrinsic nucleosome characteristics such as inclusion 
of histone variants as well as extrinsic factors such as 
ATP-dependent nucleosome remodeling, the transcriptional 
machinery and various other factors shape the dynamic 
profile of nucleosomes [24,25]. Despite the development 
of protocols and techniques that have enabled us to 
map the genome-wide nucleosome occupancy, their dy- 
namic properties are only poorly understood. Our current 
understanding of the dynamics of nucleosomes comes 
from studies performed in yeast and Drosophila. Trans- 
genic epitope-tagged histones can be inducibly expressed 
to estimate nucleosome turnover and allow detection of 
specific histone incorporation [26-28]. Alternatively, newly 
synthesized, native histones can be metabolically labeled 
with an amino acid analogue that is coupled to an affinity 
tag (CATCH-IT), which allows for detection of H3/H4 
tetramers [29]. Using these strategies, it has been shown 
that nucleosome exchange is rapid at promoters and 
coding regions, and relatively slower at heterochromatic 
regions. Furthermore, differential turnover can be quite 
localized. For instance, faster nucleosome turnover has 
been detected at Trithorax-group binding sites than at 
polycomb group protein binding sites [26,30]. 

Apart from the measurement of turnover, inducible 
expression systems with tagged histones have also con- 
tributed to our understanding of mechanistic aspects that 
pertain to histone deposition. For instance, studies from 
yeast have shown that Asfl is required for the deposition 
of H3 and that the amino termini of both H2B and H3 
are not required for their incorporation into nucleo- 
somes [27,28]. 

While these techniques measure average histone 
deposition rates across cell populations, alternative tech- 
niques such as FRAP and SNAP-tag have allowed the 
deposition of histones in individual cells to be visual- 
ized [31-33]. 

The genome-wide turnover of the histone variant H3.3 
in mammals has not been studied up until now. Global 
FRAP studies in HeLa cells with green fluorescent protein 



(GFP)-tagged histones revealed only cycling of H2B. 
Canonical H3 and H4, in contrast, exhibited very slow 
cycling, and the majority of H3 remained permanently 
bound outside S phase [33]. Slow turnover of core his- 
tones may be a feature of somatic cells since core histone 
exchange is substantially more rapid in pluripotent ESCs 
than in differentiated cell types [34]. Hence, fast turnover 
may be inherently linked to cell plasticity. 

In this study, we developed a versatile technique to 
map the dynamics of histone variant incorporation into 
chromatin in mammalian cells. Using this technique, we 
mapped the replication-independent incorporation of the 
histone variant H3.3 in mouse embryonic fibroblasts 
(MEFs). We were able to track H3.3 incorporation across 
a relatively short time window of several hours after 
induction of H3.3 as well as over a longer time frame 
of up to 72 hours. By combining our chromatin immuno- 
precipitation (ChIP) -based technique with high- throughput 
sequencing, we measured the H3.3-nucleosome turnover 
kinetics at the genome-wide level. Our results reveal three 
major categories of H3.3 nucleosome turnover: (1) rapid 
turnover at enhancers and promoters; (2) intermediate 
turnover at gene bodies; and (3) slow turnover at hetero- 
chromatic regions. We find faster H3.3 turnover at en- 
hancers and promoters is positively correlated with active 
histone modifications, including H3K4mel, H3K4me3, 
H3K9ac, H3K27ac and the histone variant H2A.Z, whereas 
slower turnover is negatively correlated with H3K27me3 
and H3K36me3 modifications. These results show that 
distinct mechanisms of histone deposition and eviction 
pertain to the dynamics of nucleosomes at different func- 
tional chromatin regions. We also show that turnover is 
related to the presence of specific histone marks, strongly 
suggesting that histone modifications are important deter- 
minants of nucleosome stability. 

Results 

An ectopic expression system to measure turnover of H3.3 

In order to track histone incorporation and thereby assay 
the genome-wide dynamics of the histone variant H3.3, we 
generated MEFs that carry a cytomegalovirus-controlled 
tetracycline (TET) transactivator and hemagglutinin (HA)/ 
FLAG-tagged H3.3 expression cassette controlled by tetra- 
cycline response elements. This TET-ON expression system 
allowed us to induce the expression of a HA/FLAG-tagged 
version of H3.3 by addition of the tetracycline analog doxy- 
cycline (DOX) (Figure 1A). In our tetracycline inducible- 
HA/FLAG-H3.3 MEF cell line, we detected robust H3.3 
expression as early as 2 hours after DOX addition that 
continued to increase until 48 hours after DOX addition. 
No tagged H3.3 expression was detected in the absence of 
DOX (Figure IB). Immunoblotting against H3.3 revealed 
that transgenic H3.3 expression levels were only a small 
fraction of those of endogenous H3.3. Furthermore we 
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Figure 1 A versatile system to study replication-independent nucleosome dynamics in mammals. (A) Schematic of T^-inducible expression 
system to study H3.3 turnover. CMV, cytomegalovirus; rtTA, reverse tetracycline-controlled transactivator; TRE, tetracycline responsive elements. 
(B) Western blot showing protein levels of transgenic HA/FLAG-H3.3 compared to endogenous H3.3. HA/FLAG-H3.3 expression 24 hours after 
DOX addition. The band marked with an asterisk is non-specific. The arrow marks transgenic HA/FLAG-H3.3. (C) Time course western blots of 
HA/FLAG-H3.3 expression. (D) Bromodeoxyuridine (BrdU) immunostaining of NIH/3 T3 cells treated with DNA polymerase inhibitor aphidicolin 
and DOX across time points of H3.3 induction. DMSO, dimethylsulfoxide. (E) Cell cycle analysis of cells treated with aphidicolin/DOX. Cells were 
stained with propidium iodide and analyzed by flow cytometry. 



verified that the HA/FLAG tags did not interfere with the 
H3K4me3 modification of H3.3 (Figure 1C). 

In order to minimize the effect of replication-coupled 
histone disassembly, we arrested the cell cycle of conflu- 
ent NIH/3 T3 MEF cells by treatment with the DNA 



polymerase (Pol) inhibitor aphidicolin. After 18 hours 
of aphidicolin treatment (time point 0) and across the 
time course of HA/FLAG-H3.3 induction, the MEF cell 
population was essentially devoid of cells in S phase 
and arrested at the Gl/S phase boundary, as indicated 
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Figure 2 (See legend on next page.) 
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(See figure on previous page.) 

Figure 2 Distribution and enrichment of H3.3 in the genome. (A) Enrichment of HA-H3.3 at various genomic regions of different annotation. 
The total read count of each category was normalized over its total length. (B) Pie chart showing the proportions and numbers of peak centers 
that fall into various genomic categories. (C) Enrichment of HA-H3.3 at various repetitive elements. AII_Rep, all repetitive elements; SAT, satellite 
regions; SINE, short interspersed nuclear element. (D) Count of histone mark peaks that overlap with HA-H3.3 peaks. Overlapped peaks were defined as 
overlapping when one peak (H3.3 or histone mark) shared 1 0% of reads with the other. ChlP-Seq libraries for various histone marks were prepared from 
cells treated with aphidicolin. (E) UCSC Genome browser view illustrating overlap of HA-H3.3 with various histone marks and H2AZ. (F) Distribution profile 
of HA-H3.3 based on mRNA expression levels as indicated by different colors. RNA-seq libraries were prepared from cells treated with aphidicolin. 
(G) Distribution profile of HA-H3.1 based on mRNA expression levels as indicated by different colors. ChlP-Seq and RNA-Seq libraries were prepared 
from untreated (no aphidicolin) cells. 



by bromodeoxyuridine and propidium iodide staining 
(Figure 1D,E). Thus, to monitor the genome-wide dy- 
namics of replication-independent H3.3 incorporation, 
we induced HA/FLAG-H3.3 expression in cells arrested 
by aphidicolin, followed by ChlP-Seq analysis using the 
HA antibody at various time points. We took a 'high- 
resolution' approach by tracking histone incorporation 
across hourly time points of early protein expression 
and across a longer time frame of up to 48 hours. 

Genome-wide characterization of H3.3 incorporation 

In order to characterize the genome-wide deposition of 
HA-H3.3, we sought to map the H3.3 distribution 72 hours 
post-induction. Consistent with previous reports from 
HeLa and mouse ESCs [1,9], we found that H3.3 is 
strongly enriched in promoters, 5' UTRs, as well as the 
3' ends of genes (Figure 2A,B,F). In terms of read density, 
promoters and TES regions were overrepresented in the 
H3.3-enriched genomic regions (Figure 2A). The majority 
of H3.3 peaks were detected at introns followed by peaks 
in intergenic regions (Figure 2B). In line with the idea that 
H3.3 carries primarily active histone modifications [13], 
we found strong overlap of H3.3 with the distribution 
profiles of histone modifications such as H3K4mel, 
H3K4me3, H3K9ac and H3K27ac. In contrast we detected 
substantially fewer regions of H3.3 incorporation at sites 
enriched with H3K27me3 (Figure 2D, 2E). Since the vast 
majority of sharp H3.3 peaks in introns and intergenic 
regions were associated with enhancer marks, including 
H3K4mel, H3K27ac and H2A.Z, we considered them as 
putative enhancers in the following analyses. 

Although most of the repetitive elements showed obvious 
depletion of H3.3, high levels of H3.3 were detected at 
telomeres, rRNA and tRNA repeats (Figure 2C). H3.3 
co-localizes with pericentromeres in several differentiated 
cell types [10]. We, however, did not find significant 
deposition of HA-H3.3 at pericentromeric regions. This 
strongly suggests that H3.3 deposition in these regions 
is a strictly replication-coupled process as we profiled 
the distribution of HA-H3.3 from cell cycle arrested and 
non-dividing cells. 

To test the relationship between H3.3 incorporation 
and gene expression, we plotted the average profiles of 
H3.3 across several groups of genes sorted according to 



their expression levels (Figure 2F). Unlike in mouse ESCs, 
we found that H3.3 is associated mostly with active genes 
and not inactive genes. 

Enrichment levels of H3.3 at TSSs, gene bodies and 
TESs were positively correlated with the transcription 
level of the associated genes (Figure 2F). Interestingly, 
we find that the relation between HA-H3.3 enrichment 
and gene expression follows a bimodal distribution at 
promoters. Only promoters of expressed genes carry 
substantial amounts of H3.3, and H3.3 enrichment levels 
do not increase further at even higher expression levels. 
In contrast, HA-H3.3 enrichment in gene body and TES 
regions was strongly positively correlated with gene ex- 
pression levels. This indicates that H3.3 in gene bodies 
may be directly transcription coupled, whereas other 
mechanisms of H3.3 deposition may exist upstream of 
TSSs. As a control alongside H3.3, we generated a MEF 
cell line that expressed canonical HA-H3.1 instead of 
HA-H3.3, which differs in five amino acids from H3.3, 
and mapped its replication-coupled deposition. Changing 
H3.3 to H3.1 abolished the defined enrichment of H3.3 
to a homogenous distribution similar to reads mapped 
from input samples, confirming its universal deposition 
throughout the mammalian genome (Figure 2E). Fur- 
thermore, the pronounced pattern of H3.3 enrichment 
around the TSS, gene body and TES was lost. Slightly 
higher levels of H3.1 were detected around the TSSs 
of inactive and low-expressed promoters, indicating 
prevalence of H3.1 at inactive genes (Figure 2G). Simi- 
lar genome-wide ChlP-Seq profiles were obtained by 
Goldberg et al. [9] with zinc finger-targeted, heterozy- 
gously tagged H3.3B in ESCs: mutating HA-H3.3 to 
HA-H3.2 or to a HA-H3.1S31 hybrid altered the H3 
distribution to a genome-wide pattern similar to that 
of input samples or ChlP-Seq profiles obtained with 
antibodies against general H3. 

H3.3 displays early and late chromatin incorporation 

The above results indicated the steady- state enrichment 
of H3.3 across the genome. In order to measure turn- 
over rates of H3.3, we profiled the dynamics of HA-H3.3 
incorporation after 0, 1, 2, 3, 4, 5, 6, 12, 18, 24 and 
48 hours of DOX treatment (Figure 3A). Using a linear 
regression model, we calculated a turnover index for 
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Figure 3 H3.3 nucleosomes display fast and slow turnover rates. (A) Genome browser view of HA-H3.3 profiles at several genomic regions 
over a 6-hour time frame. Examples of HA-H3.3 peaks that exhibit slow (left two panels), intermediate (center two panels) and fast turnover (right 
two panels) are displayed. (B) Dot plot of turnover indices calculated from two independent experiments. (C) Heat map showing relative HA-H3.3 
enrichment (red = low, green = high) across all time points sorted by turnover index from low to high. (D) Distribution of turnover indices reveals 
two populations of H3.3 nucleosome exchange. 



each individual H3.3 peak that was identified using SICER 
and found that the turnover indices are highly reprodu- 
cible from two independent experiments (Figure 3B). 

We observed very different H3.3 incorporation kinetics 
across the genome. A substantial number of sharp peaks 



appeared within 2 to 3 hours of DOX induction and 
reached their maximum enrichment well before the end 
point of the induction (Figure 3A,C). The subsequent 
decline in read numbers that was observed at peaks with 
rapid H3.3 deposition may be a result of 'read shifting'. 
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As the total read coverage becomes saturated, more 
reads may shift to peaks, which appear at later stages of 
induction. Importantly, this way, sites of enrichment 
may be identified, which would otherwise go undetected 
in a steady state. Meanwhile, we also observed many 
broad peaks, mainly in transcribed gene body regions, 
which took 12 to 24 hours to be detected (Figure 3A,C). 
Plotting the distribution of turnover rates revealed the 
existence of two populations of H3.3 peaks with slow 
turnover and fast turnover, respectively (Figure 3D). 



Distinct H3.3 exchange rates in promoters, enhancers and 
gene bodies 

To examine the H3.3 turnover rates at different genomic 
regions, we calculated the average mean and median turn- 
over indices at each region (Figure 4A). The analysis re- 
vealed that 5' UTR and promoter regions generally have 
the highest turnover indices; enhancers have the second 
highest, and 3' UTR and TES regions have the lowest 
turnover indices (Figure 4A). To examine the timing 
of H3.3 incorporation in these different regions, we 
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Figure 4 Promoters and enhancers are associated with fast H3.3 nucleosome turnover whereas gene body and transcription end site 
regions are associated with slow H3.3 nucleosome turnover. (A) Median and mean turnover indices (TIs) calculated for various genomic 
regions. Intergenic enhancers were defined as peaks that are located in intergenic regions at least 1 kb from the TSS. Intronic enhancers were 
defined as peaks located in introns only. Gene body regions were defined as broad peaks that spanned both intronic and exonic regions and 
shared significant overlap with H3K36me3. (B-E) Distribution profiles of HA-H3.3 from TSS to TES at different time points. (F) Distribution 
plot illustrating the range of turnover indices within various genomic categories. Turnover rates were log-normalized. RPKM, reads per kilobase of exon 
model per million reads. (G) Distribution plot illustrating the range of turnover indices of H3.3 nucleosomes within various genomic categories. Log of 
turnover indices was calculated and re-scaled from 0 to 1 . 
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compared the average profiles of H3.3 density of either 
all genes or active genes for the first 12 hours or after 
12 hours till 72 hours of induction (Figure 4B-E). Notable 
incorporation of HA-H3.3 was detected at the promoter 
region within 3 hours. Further increases in deposition of 
HA-H3.3 were observed within 6 hours at the TSS and 
the majority of genes reached their maximum by 12 hours 
before starting to decline (Figure 4B,D). Remarkably, 
the onset of HA-H3.3 deposition at gene bodies was 
not detected any time before the 6-hour time point 
(Figure 4C,E). H3.3 deposition at gene bodies continued 
to increase from 12 hours to 72 hours of induction, 
when it reached the steady-state level. It is interesting 
to note that the signal of HA-H3.3 at the TSS exhibited 
decreases from 12 hours to 72 hours. 

High turnover of H3.3 in promoter regions of active 
genes suggests that nucleosome turnover may be corre- 
lated with transcriptional activation. Further examination 
of turnover rates at promoters indeed showed a modest 
but positive correlation with gene expression levels 
(Figure 4F). The relationship between promoter turnover 
and gene transcription suggests additional mechanisms by 
which nucleosome exchange facilitates and/or is facilitated 
by gene activation. 

To further examine variation of the H3.3 turnover rates 
within each genomic category, we plotted the distribution 
of turnover indices separated into these categories 
(Figure 4G). The analysis revealed a relatively narrow 
range for the high turnover at promoters and 5' UTRs 
and slow turnover at gene bodies and 3' UTRs, respect- 
ively, suggesting that the H3.3 nucleosome exchange in 
these regions are controlled by distinct mechanisms at 
these respective sites (Figure 4G). Interestingly, enhancer 
regions exhibited broad variability with regards to their 
turnover rates, indicating that not all regulatory regions 
are marked by high nucleosome turnover. Instead add- 
itional factors such as histone variants or histone modi- 
fications may contribute to nucleosome stability and 
turnover. 

Rapid H3.3 nucleosome turnover is associated with active 
histone marks at promoters and enhancers 

Several studies have shown that H3.3 is enriched in tran- 
scriptional regulatory regions such as promoters and 
enhancers [1,9,15]. Consistent with these observations, 
we found that many H3.3 peaks are co-localized with 
active histone modification marks, including H3K4mel, 
H3K4me3, H3K9ac, H3K27ac and the histone variant 
H2A.Z, that are often associated with promoters and 
enhancers (Figure 2D,E) [1,35,36]. To elucidate the 
relationship between H3.3 nucleosome turnover and 
histone modifications, we sorted all H3.3 peaks based 
on their turnover rates and displayed histone modifi- 
cation signals using heatmaps. The analysis indicated 



that faster turnover rates are generally associated with 
higher levels of H3K4mel, H3K4me3, H3K9ac, H3K27ac 
and H2A.Z, whereas slower turnover is associated with 
higher levels of H3K27me3 (Figure 5A). More quantitative 
analyses confirmed that indeed H3.3 peaks with higher 
levels of active modifications are turned over faster. In 
contrast, H3.3 peaks associated with heterochromatic 
marks such as H3K9me2 and repressive marks such as 
H3K27me3 are turned over slower (Figure 5B). 

To understand the differences between the enhancers 
with fast H3.3 turnover and those with slow H3.3 turn- 
over (Figure 4G, black line), we compared the histone 
modifications associated with these two groups of enhancers 
(Additional file 1). The analysis revealed that enhancers 
that displayed high turnover carried significantly higher 
median levels of active chromatin marks whereas enhancers 
of low turnover carried higher levels of H3K27me3 
(Additional file 1). Unexpectedly, our data indicated a 
negative correlation between H3.3 turnover and H3K36me3 
levels (Figure 5A,B), in line with the idea that H3K36me3 
stabilizes nucleosomes in order to prevent cryptic tran- 
scription [37,38]. 

To evaluate the relative contribution of individual histone 
modifications to turnover rates, we compared the group of 
enhancers associated with all five active modifications to 
groups missing single histone marks. The results indicate 
that the turnover rates decreased significantly by missing 
any of the active modifications (Figure 5C). Consistent with 
this observation, association with only one of these 
modifications increased the turnover rates compared to 
the enhancers without these modifications (Figure 5D). 
In contrast to these active modifications, inclusion of 
either H3K27me3 or H3K36me3 substantially decreased 
the turnover index of H3.3 nucleosomes (Figure 5E). 

In summary, we find that H3.3 exchange rate is posi- 
tively correlated with active chromatin modifications 
(H3K4mel, H3K4me3, H3K9ac, H3K27ac and H2A.Z) 
and negatively correlated with repressive marks (H3K9me2 
and H3K27me3). While the steady-state level of H3.3 is 
positively correlated with H3K36me3, its turnover rate 
is negatively correlated with H3K36me3. 

The kinetics of HA-H3.3 incorporation at tRNA and 
rRNAs genes 

We found that HA-H3.3 is enriched in tRNAs and rRNA 
genes (Figure 2C). The incorporation of HA-H3.3 was 
apparent from 6 hours of induction and continued to 
increase till 24 hours before it started to decline at 
48 hours (Figure 6A). The incorporation of HA-H3.3 at 
rRNA genes was apparent at 12 hours and reached its 
maximum at 48 hours (Figure 6B). Both tRNA and rRNA 
exist as tandem repeat regions, are highly expressed and 
are transcribed by Pol III and Pol I, respectively. Pol Ill- 
transcribed genes share many epigenetic characteristics 
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Figure 5 Fast turnover rates are associated with active histone modifications. (A) Heat map illustrating the relation between turnover rates 
and histone marks. Turnover indices were sorted from low to high (top to bottom). Relative enrichment (red = low, green = high) levels of various 
histone marks is shown. (B) Enrichment values for various histone marks were grouped into four bins with increasing enrichment levels and 
mean turnover indices were calculated for each bin. (C) Missing any of the five 'active' modifications decreases the H3.3 turnover index of enhancer 
nucleosomes. The turnover indices for enhancer H3.3 nucleosomes containing all five active modifications or missing one single modification are 
indicated on the y-axis. (D) Association with any of the five active modifications increases the turnover index of H3.3-nucleosomes at enhancers. The 
turnover indices for enhancer H3.3 nucleosomes containing none of the five active modifications or only one modification are indicated on the y-axis. 
(E) Association with either H3K27me3 or H3K36me3 decreases the turnover index of H3.3 nucleosomes. The turnover indices are indicated on 
the y-axis for H3.3 nucleosomes containing all five active modifications (only active) or five active modifications plus either H3K27me3 or 
H3K36me3 or both. 



with Pol II transcribed genes, such as presence of 
H3K4mel/2/3 and absence of H3K27me3 [39]. Hence, 
high turnover at tRNA and rRNA repeats is consistent 
with our previous result that demonstrated high H3.3- 
nucleosome turnover at sites marked with active histone 
marks. High nucleosome exchange may directly facili- 
tate accessibility of Pol Ill-specific transcription factors 



TFIIIB and TFIIIC to provide functional promoters and 
assist in transcriptional initiation. 

Slow H3.3 replacement in telomeres 

Compared to promoters, enhancers, gene body regions, 
tRNAs and rRNAs, the incorporation of HA-H3.3 at 
telomeres exhibited the slowest rate. It was detected at 
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Figure 6 Different repetitive sequences display distinct H3.3 
turnover rates. (A-C) Mapping of HA-H3.3 to repetitive elements 
across all time points showing turnover of H3.3 at tRNA (A), rRNA (B) 
and telomeric repeats (C). 



48 hours and reached a maximum at 72 hours after 
induction (Figure 6C), suggesting that a different mech- 
anism for the HA-H3.3 deposition may take effect at 
telomeres. Heterochromatic regions are expected to have 



stable nucleosomes and most of their histone exchange is 
expected to take place during S phase. Yet, we observed 
slow rates of HA-H3.3 replacement at telomeres outside 
S phase, indicating slow replication-independent histone 
exchange that may be necessary for maintenance of 
telomeric stability. It is known that repetitive features of 
telomeres and pericentromeres have a relatively low 
affinity for nucleosomes and this may contribute to his- 
tone displacement [40]. 

Discussion 

Nucleosome dynamics and histone turnover are not well 
understood in mammals. In this report, we have devel- 
oped a robust, TET-inducible system to study the kinetics 
of genome-wide deposition of the histone variant H3.3 
in MEF cells. Using this system, we were able to measure 
turnover rates of genome-wide H3.3 target sites and were 
able to infer important differences between turnover rates 
of various genomic regions. Several categories of H3.3 
deposition kinetics were observed: invariably high turn- 
over at promoters, a broad range of turnover rates at 
enhancers, and slow turnover at gene bodies. In addition, 
repeat elements displayed dramatically different turnover 
rates with relatively high turnover at tRNA and rRNA ele- 
ments and very slow H3.3 exchange at telomeres. These 
results provide novel insights into the genome-wide turn- 
over of H3.3-containing nucleosomes and suggest distinct 
mechanisms of nucleosome assembly, stability and eviction 
in order to fulfill their function in regulating transcription 
and maintaining chromatin stability and integrity. 

Genome-wide incorporation of H3.3 

Unlike canonical histones, the histone variant H3.3 can 
be deposited onto chromatin in a replication-independent 
manner [7,8]. It is enriched in transcriptionally active re- 
gions, such as gene bodies and promoters and enhancers 
of mammalian cells [1,9]. Incorporation of H3.3 into 
chromatin destabilizes the nucleosome structure and 
may facilitate transcriptional activation by generating a 
more accessible chromatin configuration [1,2]. Its de- 
position at promoters as well as in gene bodies is asso- 
ciated with epigenetic inheritance and may contribute 
to epigenetic memories [3,16]. 

Our analysis of H3.3 distribution at the steady state 
revealed its presence at intergenic regions that overlap 
with histone modifications that typically mark enhancers, 
such has H3K4mel, H3K27ac and H2A.Z, in addition to 
the enrichment at TSS regions. In coding regions, we 
find that H3.3 is broadly distributed within the gene 
body of actively transcribed genes as well as TESs and 
regions immediately after. This association is well corre- 
lated with gene expression levels as has been reported 
from other cell lines [1,9]. At promoter regions we find 
strong differences in H3.3 enrichment between inactive 
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genes and active genes but relatively moderate differ- 
ences between different expression groups, which follows 
an almost bimodal distribution. Based on H3.3 incorpor- 
ation studies that followed H3.3 incorporation over time 
it was suggested that H3.3 incorporation at promoters 
was constitutive as only minor increases in H33 enrichment 
were observed upon interferon stimulation [16]. Our 
data argue against a direct transcription-coupled mech- 
anism of incorporation at promoters as expression levels 
and enrichment levels are bimodal, which further suggests 
that the H3.3 incorporation may function upstream of 
the transcriptional activation. Enrichment of H3.3 within 
gene body regions and at the TES showed strong correla- 
tions with transcription levels, in line with the notion that 
H3.3 incorporation in gene bodies is directly coupled to 
elongation. The ChlP-Seq profiles between gene body 
H3.3 and Pol II are highly similar. Furthermore, a physical 
interaction between HIRA and Pol II has been demon- 
strated and provides further testament for a direct link 
between H3.3 deposition and transcriptional elong- 
ation [32]. The role of H3.3 deposition in gene bodies 
remains unclear but may be implicated in conferring 
epigenetic inheritance in order to facilitate constitutive 
transcription of genes after cell replication has taken 
place. 

Dynamics of H3.3 incorporation and differential turnover 

Our results on the dynamics of H3.3 incorporation revealed 
three basic modes of H3.3 deposition kinetics. The most 
dynamic exchange of H3.3 was observed at promoters 
and enhancers. Robust signals of induced H3.3 were 
detected within 2 to 3 hours of induction at these regula- 
tory regions, which is consistent with the rapid nucleo- 
some turnover at chromatin boundaries observed in yeast 
[26]. The second category of intermediate-rate incorpor- 
ation of newly induced H3.3 was found at gene body 
regions of active genes. The incorporation of H3.3 to 
gene bodies became apparent at 12 hours and reached 
a maximum at 72 hours. The slowest incorporation of 
H3.3 was detected at telomeres, with detectable H3.3 
deposition at around 24 to 48 hours post-induction. Such 
differential turnover suggests that multiple mechanisms 
of H3.3 deposition and displacement occur, which no 
doubt will become subject to further investigation. 

The fast turnover at enhancers and promoters may in- 
volve histone chaperones such as HIRA and Atrx-Daxx 
and other chromatin remodeling enzymes as well as 
sequence-specific transcription factors. High turnover 
at promoters suggests that high nucleosome turnover is 
a pre-requisite for binding of the transcriptional ma- 
chinery as well as transcription factors. 

The concept of nucleosome versus transcription fac- 
tor competition has been well established from the 
Pho5 promoter in yeast [19,23,41]. Without stimulation, 



nucleosomes are bound to the binding sites of the trans- 
activator Pho4 and effectively repress expression of the 
gene. Only upon release of nucleosomes at the promoter 
can Pho4 bind to its target site and gene activation ef- 
fectively occur [23]. Unlike in yeast, active promoters of 
mammalian cells are occupied by nucleosomes [20]. So 
what is the role of their turnover? It is attractive to 
speculate that nucleosomes sterically compete with DNA 
binding factors and that accessibility can only be granted 
when nucleosome turnover is high [26]. 

Although simple nucleosome depletion could accom- 
plish accessibility, 'transient' nucleosome deposition may 
stabilize DNA integrity, which may be necessary for 
accurate transcription factor binding. Several families 
of transcription factors, so-called pioneer factors' [42], 
require the presence and 'guidance' of nucleosomes for 
proper binding and subsequent nucleosome remodeling 
[43]. Secondly, nucleosomes and their histone modifica- 
tions provide docking stations for a variety of transcription 
factors and chromatin remodeling factors, which may 
require continuous recruitment for initiation of transcrip- 
tion. Continuous nucleosome turnover may allow for 
accessibility and nucleosome binding for chromatin fac- 
tors at the same time, which could be facilitated by the 
presence of multiple active modifications at enhancers 
or promoters. In contrast, repression of transcription 
and epigenetic inheritance may require greater stability 
of nucleosome organization and thus slower turnover 
at these sites, which, in turn, could be facilitated by associ- 
ation with polycomb group proteins and H3K27me3 
modification [30]. 

The modest rate of H3.3 incorporation in gene bodies 
and TES regions is closely related to active transcription, 
and the accumulative signals of H3.3 in these regions 
may require multiple rounds of transcription of a gene 
unit. It is generally thought that RNA polymerase displaces 
nucleosomes during elongation, as it cannot pass by nucle- 
osomes ahead of it [44,45]. Nucleosomes must be replaced 
behind passing RNA polymerase in order to maintain the 
fidelity of transcription and avoid cryptic transcription [46]. 
Unlike other histone chaperones, HIRA possesses unique 
DNA binding properties that allow it to bind to naked 
DNA. Due to its physical association with Pol II, a 'gap 
filling' mechanism has been proposed by which HIRA 
under the guidance of PolII will fill protective H3.3- 
nucleosomes into spaces of naked DNA [32]. Our study 
reaffirms the continuous eviction and replacement of 
H3.3-containing nucleosomes during transcription. Slow 
exchange of these nucleosomes likely enhances chromatin 
stability before nucleosomes are evicted during a new 
round of transcription. H3K36me3 is placed by SET2- 
related enzymes following active transcription and plays 
a role in facilitating nucleosome deacetylation [37,47]. 
Histone deacetylation stabilizes nucleosomes and thus 
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may result in the observed negative correlation between 
H3K36me3 level and H3.3 turnover in gene body regions. 

Our analyses suggest that heterochromatic regions turn 
over more slowly than euchromatic regions. We observed 
very slow turnover at telomeres and no turnover at all 
at pericentromeric regions, both of which are known to 
be enriched in H3.3. This suggests that pericentromeric 
nucleosomes are replaced only during replication, whereas 
telomeres undergo continuous exchange of their nucleo- 
somes. Continuous exchange of nucleosomes at telomeres 
may be necessary for telomeric maintenance, although 
bulk telomeric nucleosomes may be exchanged during 
replication. 

In summary, we have provided the basis for our under- 
standing of nucleosome dynamics in mammals by pro- 
viding a detailed picture of genome-wide replacement of 
the histone variant H3.3. This study will provide the 
platform for further studies into the determinants and 
mechanisms of nucleosome exchange. 

Conclusion 

In this study, we mapped the genome-wide H3.3-specific 
nucleosome occupancy and the dynamic turnover of this 
histone variant in mammalian cells. We found that H3.3 
turnover rates vary dramatically across functionally dis- 
tinct genomic regions: turnover was highest at promoters 
and enhancers, intermediate at gene bodies and slowest 
at telomeres. Furthermore, we delineated striking correla- 
tions between turnover rates, histone modifications and 
H2A.Z, suggesting that intrinsic nucleosome properties 
such as histone modifications and histone variant inclu- 
sions are important properties of nucleosome stability. 

Materials and methods 

Cell line derivation 

cDNA of human H3.3B or mouse H3.1 (Histlh3h) was 
subcloned in-frame with a HA and FLAG tag at the 
carboxyl terminus and cloned into the lentiviral plvx- 
Tight-Puro Vector (Clontech, Mountain View, CA, USA). 
Lentiviral particles were packaged in 293 T cells with the 
psPAX2 packaging plasmid. Subsequently, we transduced 
NIH/3 T3 Tet-On @ 3G cells (Clontech) and drug-selected 
with puromycin for stable integration. 

Western blotting 

Cells were lysed with RIPA buffer, and whole- cell lysates 
were resolved on SDS-PAGE, transferred onto nitrocellu- 
lose membranes, and blotted with anti-HA (Roche, Basel, 
Switzerland), anti-FLAG (M2, Sigma, St. Louis, MO, 
USA) antibodies at 1:1,000. For anti-H3.3 western blot- 
ting, histones were isolated by acid extraction as described 
in [48] prior to immunoblotting with an anti-H3.3 anti- 
body (Abeam, ab62642, Cambridge, UK). The anti-H3.3 
antibody recognizes serine 31 of H3.3, which is not 



present on either H3.1 or H3.2, but cross-reactivity with 
other histone variants has not been tested experimen- 
tally (Abeam). 

Cell cycle analysis 

Cells were examined for BrdU incorporation with the 
BrdU Labeling and Detection kit (Roche). For further 
cell cycle and DNA content analysis, cells were fixed with 
70% ethanol, and subsequently stained with 50 (ig/ml pro- 
pidium iodide in 0.1% PBS-Triton-X with 0.2 mg/ml 
RNase A before analysis by flow cytometry. 

Time course induction and preparation of ChlP-Seq/RNA-Seq 
libraries 

NIH/3 T3 MEFs were cultured in standard conditions with 
medium containing TET-compatible fetal bovine serum. 
MEFs were grown to 100% confluence and treated with 
1 (ig/ml aphidicolin for 18 hours. Adding 3 (ig/ml doxy- 
cycline hyclate before crosslinking with formaldehyde 
at various time points induced HA/FLAG-H3.3 expres- 
sion. ChlP-Seq experiments were performed as described 
previously [49] with an antibody against HA. In order to 
correlate HA/FLAG-H3.3 turnover with the presence of 
histone modifications and/or gene expression levels, we 
prepared ChlP-Seq and RNA-Seq from cells that were 
treated with aphidicolin for 18 hours prior to crosslinking 
and cell lysis. 

Data analysis 

ChlP-Seq read mapping and peak calling 

ChlP-Seq reads were mapped to the mouse genome 
(mm9, or NCBI 37) using Bowtie [50], only allowing 
'one read per position' matching. Redundant reads were 
removed from each ChlP-Seq library. We identified ChlP- 
Seq tag-enriched peaks using SICER [51], which takes 
advantage of the enrichment information from neigh- 
boring bins to identify spatial clusters of signals that are 
unlikely to appear by chance. In this way, we could iden- 
tify the dispersed epigenomic domains such as broad 
H3.3 domains on gene bodies. For calling of H3.3 peaks, 
we set the window size to 200 bp and gap size to 600 bp. 
Peaks with a false discovery rate higher than 0.01 were 
excluded from the analysis, compared with the input 
control libraries. 

Analysis of H3.3 enrichment in genie regions and repetitive 
elements 

The RefSeq gene annotation (mm9 or NCBI37) was down- 
loaded from the UCSC browser. The mouse genome 
was classified into 5' UTR, exon, intron, 3' UTR and 
intergenic based on the UCSC annotation. Then we de- 
fined another two categories: promoter (3 kb upstream 
of the TSS) and TES (from TES to 5 kb downstream). 
The total mapped reads in each category were normalized 
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to its total length. All mouse consensus repetitive element 
sequences, including pericentromeric element, LI, SINE1 
and telomeric element, were obtained from Repbase Up- 
date (RU) [52]. The relative enrichment of H3.3 in each 
category of repetitive element was calculated by the pro- 
portions of reads mapped to repetitive elements normal- 
ized over input. 

H3.3 peak annotation and distribution profiling 

We used the position of the center of each peak for peak 
annotation. Two peaks from different libraries were con- 
sidered as overlapping if the overlapped region accounts 
for 10% of the length of one peak. We classified genes 
with RPKM (reads per kilobase of exon model per million 
reads) <1 as lowest (inactive) genes and the other genes 
as active genes. We grouped genes into five categories 
with equal numbers of genes (highest, high, intermediate, 
low, lowest) based on their RPKM values. For each gene, 
reads falling in peaks were counted according to their 
shifted positions in 50 bp windows for the regions from 
3 kb upstream of the TSS to the TSS and from the TES 
to 5 kb downstream of the TES. Within gene bodies, 
reads falling in peaks were counted according to their 
shifted positions in windows equal to 1% length of each 
gene. The number of reads in each window was normal- 
ized by the total number of bases within the window, and 
the total number of peak-filtered reads in the correspond- 
ing sample to obtain a normalized read-tag density. 

RNA-Seq data analysis 

The reads from RNA-Seq libraries were mapped to the 
mouse genome (mm9, or NCBI37) using TopHat [53], a 
fast splice junction mapper. The gene expression level 
was measured by RPKM. 

Quantifying peak dynamics 

The total number of H3.3 peaks was identified by peak- 
calling at the 72-hour time point when both early appear- 
ing and late appearing peaks were readily detected. The 
number of reads in each of these peak regions was re- 
corded and normalized over the total mapped reads for 
each ChlP-Seq library. The relative H3.3 enrichment of 
each peak was calculated by normalizing the normalized 
reads in the peak over the normalized reads in input. A 
linear regression model was used to calculate turnover 
indices (77) for each individual peak. Assuming that en- 
richment of H3.3 at 0 hour is E 0 , then E t = 27 x t + E 0 , 
where E t equals H3.3 enrichment at each time point (E lf 
E 2 , E 3 ,...E n ), t = time point (1, 2, 3,...n). For peaks that 
reached their maximum enrichment before the end time 
point of analysis (t = 72 h), multiple linear regression co- 
efficients were calculated by fitting the end time points 
from time point of maximum enrichment to 72 h and 
t was adjusted correspondingly. We adopted the regression 



coefficient with the best fit (smallest P-value) as the 
turnover index of the peaks. The turnover index was 
scaled from 0 and 1 in order to compare the reproduci- 
bility between duplicate experiments. 

Accession numbers 

Our ChlP-seq and RNA-seq data sets have been deposited 
in the Gene Expression Omnibus data base with accession 
number GSE51505. 
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Additional file 1: Figure SI. Fast turnover of enhancer H3.3 
nucleosomes are associated with higher active histone modification and 
lower repressive modification. The H3.3-containing nucleosomes at 
enhancers were separated into two groups (fast turnover and slow 
turnover) based on the distribution of turnover indices in Figure 4G. The 
histone modification levels of the two groups are displayed and 
compared using box plots. 
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